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Abstract 

Convergence properties of binary stationary subdivision schemes for curves have been analyzed using 
the techniques of z-transforms and eigenanalysis. Eigenanalysis provides a way to determine derivative 
continuity at specific points based on the eigenvalues of a finite matrix. None of the well-known subdivi- 
sion schemes for curves have complex eigenvalues. We prove when a convergent scheme with palindromic 
mask can have complex eigenvalues and that a lower limit for the size of the mask exists in this case. 
We find a scheme with complex eigenvalues achieving this lower bound. Furthermore we investigate this 
scheme numerically and explain from a geometric viewpoint why such a scheme has not yet been used 
in computer-aided geometric design. 
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1 Introduction and Notation 

Subdivision schemes have been used successfully in computer-aided geometric design for more than two 
decades. The basic idea is to consider a mesh of points in either one or two dimensions and to construct a 
curve or a surface out of this mesh by an iterative procedure. This allows the representation of geometric 
objects by a relatively small number of points. In this paper we only consider subdivision schemes for curves. 
In the following paragraph, we review the required background very briefly and refer to [5] and [6] for a more 
detailed introduction. 

We start with an inital mesh of points = and associate with it a sequence of control points 

P*^ = {P^}. As an example consider the case k — 1; then we can view the previous construction as a 
bi-infinite sequence of points 

{. . . , (-1, Pi,), (-1/2, Pi,), (0, Po'), (1/2, Pt), . . .} 

A binary subdivision scheme doubles the number of points at each step, obtaining all new points only 
using the points of the last step. The mesh Nk is refined to Nk+i for primal schemes and to a mesh 
|2-(fe-i)^ - 2-('=-2)^} foj. ^^ai schemes (see Figure [IJ. 



To simplify notation we simply suppress this underlying refinement of the mesh and focus on the con- 
struction of the sequence of control points P''. The main data defining a binary stationary subdivision 
scheme is given by a so-called mask a — {ai G M.\i G cr(a) C Z} and refinement rules 
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where (T(a) is called the support of Sg,. The iteration of a subdivision scheme is denoted by p'^+^ = SaP''. 
Furthermore we consider only palindromic masks with finite width, i.e. we suppose (j{sl) — {—n,—n + 
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Figure 1: Difference between primal and dual refinement of the mesh 

1, • • • ,n — l,n} and = a_i if n is odd (primal schemes) and a{a) — {— n + 1, • • • ,n — l,n} and = ai-t 
if n is even (dual schemes). 

The next natural step is to analyze the limiting properties of a scheme as fc — > cx). Let p''{t) denote the 
function obtained from the control points P*^ (attached to the underlying mesh) by linear interpolation. A 
subdivision scheme Sg^ is termed C™-convergent if there exists a C™ function / : M — > R such that 

hm max ||/(t)-p'=(i)|| -0 

fe— >oo A<t<B 

for every fixed A,BeM. and if f{t) ^ for some initial sequence of control points. Two paradigms of 
analysis have been developed to answer whether a scheme is convergent. First, we define the z-transform 
of a subdivision scheme by 

Se. ^ Sa(^)- (2) 
'ieo-(a) 

giving a one-to-one correspondence between the mask a and a formal Laurent series. We have the following 
well-known results (see [1] or [2]): 

Theorem 1.1. Sa is not convergent if Sa{l) ^ 2 or So(— 1) ^ 0. 

Hence convergent schemes have Sa(z) = (1 -I- z)sb{z) where the scheme S'b relates to via A(S'aP*') = 
S'bAP''^ where AP'' = - Pt~'^\i € Z}. We call S'b the difference scheme associated to S^. 

Theorem 1.2. Sa is C'^ -convergent if Sb is contractive, i.e. Sb maps any initial sequence of control points 
to upon iteration. 

Note that contractivity of S'b can be checked by computing if ||Sb|| < 1 where 

ll^bll =maxJ^|62,|,^|62,+i| 



This reduces the question of convergence to a direct computation. Theorem 11.21 can be improved to an "if 
and only if" statement in the case of asymptotic equivalence of the operators S'a and Sb (see ^ for details). 
The easiest way to improve the level of differentiability is given in the next theorem. 

Theorem 1.3. Let Sa{z) = ^^-^^Sg{z). If Sq is C'^ -convergent then Sa is C'^'^^ -convergent. 

Hence we can try to divide out -^^^^ factors from the z-transform and check the resulting scheme for 
convergence. Conversely we can improve the level of smoothness of a convergent scheme by multiplying 
the z-transform by ^^^-^^ factors. Note that the correspondence ^ provides the corresponding mask a. 
Expositions and proofs of the previous results are found in ^2, and |3j. 
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2 Complex Eigenvalues: An Example and Minimum Width 



An alternative way of analyzing convergence is provided by setting up a matrix A associated to Sa_. Note 
that the refinement rules ([T]) require an infinite matrix but to investigate smoothness near a given point 
it suffices to perform eigenanalysis on a finite matrix A called the local subdivision matrix. Eigenanalysis 
provides an upper bound on the level of derivative continuity. If the width of 5a is n then A e R"^". For a 
more detailed dsicussion of eigenanalysis consider [3] or [S]. 

We ask the question if the matrix A can have complex eigenvalues for a convergent scheme S^- Reviewing 
many well-known subdivision schemes for curves the surprising result is that all classical schemes such as B- 
spline schemes or the four-point scheme (see ^3^ or [5^) do not have complex eigenvalues. A priori there is no 
reason for a matrix with real entries not to have complex eigenvalues. We remark that complex eigenvalues 
have been observed in the analysis of some subdivision schemes for surfaces such as the V3-scheme (see [1]). 
For curves no complex eigenvalue scheme seems to be known. Therefore we ask three questions: 

1 . Does there exist a convergent binary subdivision scheme with palindromic mask having complex eigen- 
values? (Answer: Yes!) 

2. What is the minimum width of such a scheme? 

3. Why did complex eigenvalue schemes not appear "naturally" in computer-aided geometric design for 
curves? 

Proposition 2.1. There does not exist a binary subdivision scheme Sa with width less than or equal to 5 
and complex eigenvalues in the local subdivision matrix A. 

Proof. We begin with the case n — 5. First note that since Sa(z) = a2Z^ + aiz + ao + aiz^^ + a2Z^^ we 
get using Theorem 11.11 that 2 — Sa(l) = 2a2 + 2ai + oq and = Sa(— 1) = 2a2 — 2ai + ao for a convergent 
scheme. This implies the conditions: 

oq = 1 — 2a2 and ^ 2 
Therefore the local subdivision matrix A is given by 
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which has eigenvalues {1, ^, ^ ~ 2a, a, a} and since a G M we conclude that no convergent scheme with a 
mask of width 5 can have complex eigenvalues. The cases n < 4 are dealt with in the same way and all 
matrices have real eigenvalues in these cases. □ 

Note that the previous analysis fails for the case n = 6. We show now that it is possible to find a 
convergent scheme of width 6 with "simple" coefficients having a small common denominator. 

Proposition 2.2. The mask a = {a_2 = —jQ,a-i = ^,ao = fi'^^i = t''^2 — '^1^2 = ^Jo^ defines a 
-convergent binary scheme Sa with complex eigenvalues in the local subdivision matrix A. 
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Proof. Since s^{z) = (1 + z) (-f^ + ^ + | + A _ (i + z)sb{z) we find |15b(z)|l = | < 1 so that 

Sn is C'^-convergent by Theorem 11.21 The local subdivision matrix is given by 



which has eigenvalues {- 
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□ 



This means that complex eigenvalues impose conditions on the width of a mask. It is in general desirable 
to have masks with relatively small size to guarantee locality, i.e. moving a control point does not affect 
parts of the curve far away from this point. Many classical schemes have width between 6 and 10. Therefore 
we can conclude that there is no algebraic reason why we have not yet seen a complex eigenvalue scheme in 
practice. 

We remark that we cannot use Theorem 11.31 to find a C^-convergent scheme with width 6. Since 
s^{z) = (2az^ + 2(5 -a)z + 2-Ab+ + ft ) = ^^^Sq{z) using Theorem O we get the dif- 

ference scheme Sq. Using Theorem II. II on Sq gives 



= sq(-l) = 2(1 + 4a -46) 



This condition implies that the eigenvalues of the local subdivision matrix A are real. But obviously we can 
use Theorem 11.11 to construct a scheme with higher smoothness out of the complex eigenvalue scheme we 
presented by multiplying the z-transform with -^^-^ factors. For example we get from the scheme given in 
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It again has a pair of complex conjugate eigenvalues. 

3 Numerical Comparison: Complex vs. Real Eigenvalues 

Our goal is to answer the last of our three questions and we start by comparing the results of several schemes: 





Mask 


Smoothness 


(a) 


{a-2 = 


-1/10, a- 


-1 = 3/10, ao = 4/5, ai = 4/5, as = 3/10, aa = -1/10} 






(b) 


{a~3 


-1/20, a. 


-2 = 1/10, a_i = 11/20, ao = 4/5, ai = 11/20, aa = 1/10, ag = 


-1/20} 




(c) 


{fl-i = 


1/2, ao = 


l,ai = 1/2} 




c° 


(d) 


{a_2 = 


1/8, a_i = 


= 4/8, ao = 6/8, ai = 4/8, aa = 1/8} 







Scheme (d) is derived from cubic B-splines and scheme (c) is known as "simplest" or "two-point" scheme; 
it is easy to check that (c)-(d) have only real eigenvalues. We compare the schemes with the most important 
test sequence of control points given by 



{0,0,0,0,1,0,0,0,0} 
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on the integer mesh iV" on [—4,4]. The reason for its importance is that ah convergence properties can be 
derived by investigating this sequence and using hnearity of the scheme (see [3]). The schemes were iterated 
10 steps and hnear interpolation was used to display the function constructed from P^". The results are 
shown in Figure [21 




Figure 2: (a)-(b) complex eigenvalue schemes; (c)-(d) two "standard" schemes. 

The cubic B-spline scheme is and produces the expected smooth curve. The simplest scheme produces 
the classical C" tent function. For the complex eigenvalue scheme (a) we find a C° curve, which is neither 
smooth nor contains clear connecting points between several parts of the curve like displayed for the tent 
function. Instead we see small oscillations or "ripples". The complex eigenvalue scheme (b) gives a C"^ 
curve with two small negative bumps outside [—2, 2]. In computer-aided geometric design undesired ripples 
or bumps of the limiting curve are usually called artefacts. In the last two decades the focus has been to 
find schemes which produce smooth curves or curves with well-defined sharp transitions without artefacts. 
This is one possible explanation why practioners might have overlooked or dismissed schemes with complex 
eigenvalues. In the next section we show why artefacts are expected to appear. 

4 A Geometric Explanation using Dynamical Systems 

We aim to explain why we expect artefacts in every subdivision scheme for curves with complex eigenvalues. 
Fix a convergent scheme Sa with one pair of complex conjugate eigenvalues. Fix a mesh point Xk^ G . 
Consider the n control points in P*^" attached to mesh points closest to Xk^. By shifting the indices we can 
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assume without loss of generality that these points are given by Pj " for i S {1,2, ... ,n}. Define a vector 

z;fc„ = (Pf°,P^,...,P^) (3) 

This process can be continued for the n-points closest to Xk^ in p'^o+^, pfeo+2^ ^^^^ Note that the local 
subdivision matrix A associated to 5a produces this sequence of vectors under iteration: 

Vk+i = Avk for k > kg and Vk^ given by ^ 

This means we are dealing with a particularly simple discrete dynamical system defined by the map A : 
M" M" . It is known from eigenanalysis that the eigenvalues fii of A must obey the restriction fii — 1 and 
< 1 for i e {2, 3, . . . , n} (see [S]). Furthermore it is immediate by the convergence of that 

Avk — Vk+i — > w as fc ^ oo for some vector vector v with identical components. 

In fact we find that Av — v so that w is a fixed point. Since we have at least one pair of complex conjugate 
eigenvalues there exists a change of basis so that A has a two-dimensional eigenspace R in which A is a 
rotation-scaling with scaling factor A < 1. The scaling factor is given by the real part of the pair of complex 
conjugate eigenvalues. Due to the rotational component the convergence to v is not monotone. This lack 
of monotonicity is reflected in the different behaviour of the complex eigenvalue schemes in comparison to 
some of the classical schemes. 

In our examples of complex eigenvalues there is another reason for non- monotone convergence caused by 
the two real negative eigenvalues (see Proposition l2 . 2p . This gives alternating convergence of two components 
of Vk to V. This alternating convergence is unavoidable for complex eigenvalue schemes with width 6 as the 
next proposition shows. 

Proposition 4.1. Every binary complex eigenvalue subvision scheme with palindromic mask of width 6 has 
at least 2 eigenvalues with negative real part. 

Proof. The local subdivision matrix A is given by 

6 \ 

c a 
c 6 
6 c a 
a c 6 
Q h c a ) 

Using the condition a + fe + c= lwe first eliminate c and then compute the eigenvalues explicitly. They are 
given by: 

/ii = 1, /i2,3 = a, /i4 = 6 — a, /i5^6 = 1 — a — 6 ± \/l + 2a— la^ — 66 + 2a6 + 96^ 

For the eigenvalues to be complex we get that the terms inside the square-root for /15.6 must be negative. 
The inequality to be satisfied is: 

1 + 6-2^2(1-56+862) 
a< ^ 

But we can calculate directly that 

1 + 6-2^/2(1 - 56 + 862) <0 for 6 e K 
Therefore we can conclude that /i2,3 = a < 0. □ 

We leave it as an open question whether any binary complex eigenvalue scheme with palindromic mask 
must have eigenvalues with negative real part. Hence it might be possible to avoid alternating convergence 
for a complex eigenvalue scheme but non-monotonicity can obviously not be avoided. Hence we expect 
artefacts to appear. 
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5 Conclusion 



We have demonstrated that subdivision schemes for curves can have complex eigenvalues. We have shown 
that the minimum width of a scheme with complex eigenvalues is 6. Furthermore we defined a scheme 
achieving this minimum with simple coefficients. Numerical comparison of complex eigenvalue schemes with 
two well-known schemes showed that the main reason why complex eigenvalue schemes have been overlooked 
or disregarded is the appearance of ripples or artefacts. They appear due to the lack of monotonicity of the 
dynamical system defined by iteration of the local subdivision matrix A. 

We recommend to consider the use of complex eigenvalue schemes in computer-aided geometric design if 
there is a need for curves or objects with rough surfaces. If the scaling factors of complex conjugate eigenval- 
ues and the magnitude of real negative eigenvalues is small we expect a limiting curve with small artefacts. 
This can be used to scale how pronounced the artefacts should appear. Note that we can always construct a 
surface subdivision scheme from a complex eigenvalue curve scheme by using the tensor product of the mask. 
Therefore we can hope to use these schemes to define geometric objects which do not appear completely 
smooth or with sharp transitions and close a gap in the current practical applications of subdivision schemes. 
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